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ABSTRACT 

We present a simple metallicity estimator based on the logarithmic [Nil] A6584A/Ha 
ratio, hereafter N2, which we envisage will become very useful for ranking galaxies 
in a metallicity sequence from redshift survey quality data even for moderately low 
spectral resolution. 

We have calibrated the N2 estimator using a compilation of H ii galaxies having 
accurate oxygen abundances, plus photoionization models covering a wide range of 
abundances. The comparison of models and observations indicates that both primary 
and secondary nitrogen are important for the relevant range of metallicities. 

The N2 estimator follows a linear relation with log(0/H) that holds for the whole 
abundance range covered by the sample, from about l/50th to twice the Solar value 
(7.2 < 12+log(0/H) < 9.1). We suggest that the ([Sii] AA6717,673lA/Ha) ratio (here- 
after S2) can also be used as a rough metallicity indicator. Because of its large scatter 
the S2 estimator will be useful only in systems with very low metallicity, where [N ii] 
A6584A is not detected or in low resolution spectra where [Nil] A6584A is blended 
with Ha. 

Key viTords: galaxies: abundances - galaxies: stellar content - galaxies: evolution 



1 INTRODUCTION 

, The key for accurate determination of oxygen abundance 
' in gaseous ionized nebulae is a precise measurement of the 
weak auroral forbidden emission line [O III] A4363A. In nar- 
, row emission line starforming galaxies the temperature sen- 
sitive [O III] A4363A line intensity correlates with the overall 
abundance, being relatively strong in very low metallicity 
systems (12-|-log(0/H) < 7.8) and becoming undetectable 
even for moderately low metallicity galaxies (12-|-log(0/H) 
> 8.3). As a result, for most of the starforming galaxies 
[O III] A4363A is unmeasurably weak. 

For the large majority of starforming regions the oxygen 
abundance is therefore estimated using empirical methods 
based on the relative intensities of strong, easily observable, 
optical lines. Although abundances derived in this way are 
recognized to suffer considerable uncertainties, still they are 
believed to be able to roughly trace general trends in galax- 
ies. The most widely used empirical abundance calibrators 
are the R23 (Pagel et al. 1979) and, recently, the 823(4) 
(Vi'lchez & Esteban 1996; Di'az & Perez-Montero 2000; Oey 
& Shields 2000) parameters. 

* Visiting Professor, INAOE, Puebla, Mexico 
t Visiting Fellow, loA, Cambridge 



The R23 method was first proposed by Pagel et al. 
(1979) and subsequently developed and calibrated by many 
authors (Edmunds & Pagel 1984, McCall, Rybski & Shields 
1985, Dopita & Evans 1986, Torres-Peimbert, Peimbert & 
Fierro 1989, McGaugh 1991). It is defined as the sum of the 
flux of [O 11] A3727A and [O iii] AA4959,5007A Unes relative 
to H/3 (R23 = ([On] A3727A + [Oiii] AA4959,5007A)/H/3). 
The R23 method can be used to estimate abundances up 
to relatively high redshifts, but there are a few problems 
associated with the use of this estimator. Firstly, it is bi- 
valued, i.e. a single value of R23 can be due to two very 
different oxygen abundances. Secondly, a very large fraction 
of the starforming regions lie on the ill-deflned turning zone 
around 12-|-log(0/H) ~ 8.1 where regions with the same R23 
value have oxygen abundances which differ by almost an or- 
der of magnitude. Thirdly, the R23 method requires spec- 
trophotometric data and given the wavelength range cov- 
ered, the reddening correction of the lines involved becomes 
crucial. Finally a characteristic which is readily apparent 
is the large scatter in the R23 vs. oxygen abundance cali- 
bration, larger than accounted for by observational errors 
(Kobulnicky, Kennicutt & Pizagno 1999). 

The S23 parameter introduced by Vflchez & Este- 
ban (1996) is deflned as the sum of the flux of [Sii] 
AA6717,673lA and [S lll] AA9069,9532A relative to H/3 (S23 
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= ([Sii] AA6717,673lA + [Siii] AA9069,9532A)/H/3). Un- 
like the R23 method, the relation between S23 and oxygen 
abundance remains single valued up to a metallicity slightly 
higher than solar. Diaz & Perez-Montero (2000, hereafter 
DPOO) show an empirical calibration of the S23 parameter 
with a somewhat reduced scatter as compared to that with 
R23. This calibration should be improved by the inclusion 
of high quality data at both the low and the high metallic- 
ity ends, but otherwise looks very promising for metallicities 
up to solar. On the other hand, Oey & Shields (2000) ar- 
gued that the S23 method is more sensitive to the ionization 
parameter than R23, and propose to include the emission 
of [Siv] to overcome the limitations of S23, introducing the 
S234 method. A problem with this estimator is that good 
data is still scarce, and that the detection of the sulphur 
lines in the near-infrared is limited to galaxies with redshifts 
smaller than about 0.1. 

In all, the estimation of abundances of large number 
of galaxies, particularly in the slightly subsolar to oversolar 
range, is still a difficult problem. This metallicity regime 
contains most of the H 11 regions in early spiral galaxies and 
the inner regions of most late type galaxies, therefore the 
description of the metallicity distribution in galaxies cannot 
be complete without it. 

Storchi-Bergmann et al. (1994) suggested the use of N2 
= [Nil] A6584A/Ha , as an abundance estimator. Their 
calibration of the N2 vs. 0/H relation was improved by 
Raimann et al. (2000) that proposed a new calibration using 
the Terlevich et al. (1991, hereafter T91) sample. 

In the present work we have used the best available 
abundance data for starforming galaxies in combination 
with photoionization models to explore the usefulness of 
the N2 abundance estimator and to calibrate it in terms 
of metallicity. 

The data set is presented in Section 2. The new cali- 
bration of the N2 estimator with photoionization models is 
discussed in Sections 3 and 4. The concluding remarks are 
given in Section 5. 



2 THE DATA 

We have compiled from the literature a sample with the best 
available data on the lines of interest. 

The published data included in our sample was selected 
for having high signal-to-noise, spectral resolution better 
than ~ 8 A in order to separate [Nil] A6584A from Ha, eas- 
ily accessible line strengths and errors. In computing abun- 
dances and line ratios, we have adopted the c(H/3) reddening 
factor published by the respective authors. 

The sample of low metallicity galaxies is based on a 
compilation by Kobulnicky & Skillman (1996) and is not 
intended to be complete, but includes the majority of the 
best spectrophotometric data available for H 11 galaxies. 

The data for the metal-poor galaxy sample (i.e., with 
available [O ill] A4363A intensity and 12+log(0/H) < 8.4) 
are presented in Table 1, together with their correspond- 
ing observational errors propagated from the emission line 
uncertainties quoted in the original references. We re- 
determined the electron temperatures and abundances ap- 
plying the method and expressions from Section 4 of Pagel et 
al. (1992) to the reddening corrected emission line fluxes. In 



Table 1, objects within references 1-17 had their metallicities 
computed using the [O lll] A4363A line intensity. Galaxies 
with more than one H 11 region may appear more than once 
in the table. The comparison of the tabulated values can 
give an idea of observational and systematic errors and pos- 
sible spatial variations. We have added to the sample from 
the literature, our measurements of six metal-poor galaxies 
discovered by us on the Anglo- Australian 4.0m telescope in 
August 1996 and August 1997 (Terlevich et al. 2001, here- 
after Paper II). The total number of low metallicity objects 
with 12-|-log(0/H) computed using the [Olll] A4363A line 
intensity is 108. 

The sample of metal-rich galaxies (128) includes data 
from DPOO; Castellanos, Di'az & Terlevich (2001, hereafter 
CDTOl) and T91. For galaxies with weak or none [O ill] 
A4363A emission (i.e., 12-flog(0/H) > 8.4), we derived the 
metallicity using R23 (the upper branch analytic expression 
by McGaugh, published in Kobulnicky et al. 1999) or S23 
(Di'az & Perez-Montero 2000) methods, depending on the 
availability of the [S ill] AA9069,9532A lines. The errors for 
the metallicity values derived by the R23 parameter already 
account for the intrinsic 0.2 dex uncertainty of the method. 
The 55 metal-rich galaxies included in Table 1 with refer- 
ences 18-26 correspond to the compilation in DPOO and to 
CDTOl and had their metallicities estimated by the R23 or 
S23 methods. The rest of the metal rich galaxies, as they are 
from T91, were not included in the table. 

The whole range of oxygen abundances covers from ~ 
2% solar (Solar is taken as 12-|-log(O/H)0 = 8.91) for IZwl8 
to more than solar for some regions for which detailed mod- 
eling has been performed (see compilation in DPOO). 



3 RESULTS 

The N2 parameter is defined as 

N2 = log([Nll] A6584/Ha). (1) 

The relation between N2 and the oxygen abundance 
is shown in Figure 1. It can be seen that these two pa- 
rameters are well correlated (linear correlation coefficient 
of 0.85) and that a single slope is capable of describing the 
whole metallicity range, from the most metal-poor to the 
most metal-rich galaxies in the sample. The heavy solid line 
represents the linear fit to the N2 vs. 12-|-log(0/H) relation. 
Least squares fits to the data simultaneously minimizing the 
errors in both axes, give 

12+log(0/H) = 9.12(± 0.05) + 0.73(± 0.10)-N2 , (2) 

The thin solid line in Figure 1 shows the linear fit from 
Raimann et al. (2000). Given the differences among our sam- 
ples and the smaller number of points in Raimann et al. 
data, we conclude that both solutions are in basic agree- 
ment. However both fits dramatically disagree with the cal- 
ibration proposed by Storchi-Bergmann et al. (1994, section 
4). 

3.1 Photoionization models 

We have computed a grid of photoionization models using 
CLOUDY (Ferland 1996) for abundances [0/H] = 0.04, 0.2, 
1 and 2. All the elements apart from N are assumed to be 
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Table 1. Chemical abundances and [Nil] /Ho ratios. 



Galaxy 


12+log(0/H) 




N2 




Ref. 


Galaxy 


12+log(0/H) 




N2 




Ref. 


I Zw 18nw 


7.19d±0 


046 


-2 


527+0 


068 


1 


1437+370 


7 


971+0 


064 


-1 


678+0 


022 


4 


I Zw 18se 


7.309±0 


057 


-2 


317+0 


017 


1 


Pox 139 


7 


983+0 


083 


-1 


669+0 


123 


12 


SBS 0335-052W 


7.285±0 


051 


-2 


100+0 


119 


2 


UM 462 


7 


984+0 


063 


-1 


676+0 


059 


10 


UGCA 292-1 


7.320±U 


050 


-1 


943+0 


046 


3 


N5253 A 


7 


991+0 


069 


-1 


106+0 


010 


10 


UGCA 292-2 


T or* "1 In 

7.361±U 


060 


-2 


104+0 


060 


3 


N5253-6 


8 


025+0 


066 


-1 


053+0 


020 


13 


0940+544 N 


7.o77±U 


040 


-2 


224+0 


055 


4 


N5253-5 


8 


094+0 


067 


-1 


123+0 


021 


13 


HS 0822-1-3542 




045 


-2 


260+0 


040 


5 


N5253-1 


8 


144+0 


076 


-1 


182+0 


012 


13 


1159+545 


7.473±0 


045 


-2 


067+0 


040 


4 


N5253-2 


8 


167+0 


072 


-1 


215+0 


021 


13 


1415+437 


T con 1 n 


047 


-1 


922+0 


014 


6 


N5253-4 


8 


180+0 


069 


-1 


179+0 


024 


13 


0832+699 


V rr o 1 n 


052 


-1 


982+0 


015 


4 


N5253 B 


8 


266+0 


104 


-1 


161+0 


036 


10 


T1214-277 


T crn^ 1 n 

7.59d±U 


052 


-2 


493+0 


048 


7 


T1304-386 


8 


001+0 


067 


-1 


508+0 


031 


10 


UGC4483 S 


'7 cciA In 

7.d01±0 


055 


-1 


982+0 


019 


8 


UM 469 


8 


001+0 


080 


-1 


181+0 


045 


10 


1211+540 


^7 j::? o T 1 n 

7.d87±U 


054 


-2 


142+0 


022 


4 


37-27 


8 


012+0 


064 


-1 


630+0 


025 


9 


1249+493 


"7 '70 1 1 n 

7.7zl±U 


053 


-1 


926+0 


054 


6 


1533+469 


8 


012+0 


062 


-1 


261+0 


Oil 


6 


537-69 


V rr 1 n 


101 


-1 


628+0 


046 


9 


Mrk 600 


8 


019+0 


068 


-1 


873+0 


122 


7 


T1304-353 


^ ^ /I r" In 

7.745+0 


055 


-2 


263+0 


116 


10 


1135+581 


8 


019+0 


069 


-1 


625+0 


006 


4 


70-05b 


7.751+0 


087 


-1 


719+0 


086 


9 


Pox 108 


8 


026+0 


089 


-1 


669+0 


123 


12 


UM 461 


T '7n c 1 n 

7.795±0 


057 


-2 


251+0 


022 


7 


Pox 4 NW 


8 


031+0 


082 


-1 


675+0 


123 


12 


C1543+091 


7.79d±0 


059 


-1 


939+0 


109 


10 


0946+558 


8 


034+0 


067 


-1 


645+0 


014 


4 


1331+493 N 


T o 1 n 1 n 

7.8l9±0 


057 


-1 


940+0 


027 


4 


0948+532 


8 


039+0 


068 


-1 


617+0 


020 


4 


1331+493 S 


7.911+0 


062 


-1 


441+0 


031 


6 


T1345-420 


8 


049+0 


067 


-1 


716+0 


048 


10 


1152+579 


T o r* n 1 n 

7.850+0 


057 


-1 


844+0 


022 


4 


II Zw 70 


8 


067+0 


090 


-1 


346+0 


029 


11 


Mrk 36 


7.865+0 


099 


-1 


827+0 


087 


11 


T1334-326 


8 


088+0 


072 


-1 


928+0 


106 


10 


25-10 


7.865+0 


076 


-1 


140+0 


154 


9 


II Zw 40 


8 


104+0 


077 


-1 


708+0 


024 


10 




7.868+0 


082 






123 


12 


TT 7^f.r AC\ 

il ZjW 


8 


183+0 


091 






UUJ- 


Y\ 


Mrk 36 


7.872+0 


064 


-1 


901+0 


100 


10 


N6822-HuX 


8 


104+0 


086 


-1 


740+0 


100 


14 


46-17 


7.891+0 


061 


-1 


968+0 


028 


9 


T1004-296 SE 


8 


117+0 


072 


-1 


290+0 


020 


10 


C1148-203 


7.896+0 


063 


-1 


775+0 


029 


10 


T1004-296 NW 


8 


204+0 


083 


-1 


386+0 


016 


10 


Pox 105 


7.904+0 


081 


-1 


669+0 


123 


12 


Tol 35 


8 


121+0 


087 


-1 


547+0 


123 


12 


297-24 


7.939+0 


081 


-1 


325+0 


083 


9 


T0633-415 


8 


144+0 


072 


-1 


373+0 


150 


10 


C0840+120 


7.940+0 


063 


-1 


615+0 


056 


10 


T0633-415 


8 


144+0 


071 


-1 


391+0 


008 


7 


Fairall 30 


7.949+0 


064 


-1 


634+0 


022 


10 


Fairall 2 


8 


169+0 


074 


-1 


155+0 


066 


10 


Tol 2 


7.965+0 


080 


-1 


351+0 


123 


12 


T1324-276 


8 


170+0 


077 


-1 


895+0 


024 


10 
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1994; (9)Terlevich et al. 2001; (lO)Campbell, Terlevich & Melnick 1986; (ll)Garnett 1990; (12)Kunth & 
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primary and are therefore scaled as oxygen. For nitrogen we 
used [N/O] = 0.08 + [O/H] that represents a realistic combi- 
nation of early primary N plus later secondary N and repro- 
duces the observed behaviour of N/O with O/H (Henry, Ed- 
munds & Koppen, 2000 and references therein). The square 
brackets in the expressions denote the value relative to that 
of Orion, i.e., N: 7.00 x 10"^ O: 4.00 x 10"". For simplic- 
ity we have used single star models with Te// = 45000 K. 
Changing to cluster models or including a metallicity depen- 
dent Te/ / do not change the basic results. The photoioniza- 
tion models are shown as dashed and dot-dashed lines in 
Figure 1. The lines join models with the same ionization 
parameter ( U) . The whole range of the data is comprised 
between the models with log (U)=-1.5 and log (U)=-3.0. It 
can be seen that within the errors there is a good superpo- 
sition between models and observations. 

Models with either pure primary or pure secondary ni- 
trogen fail to reproduce the slope of the relation. The change 
in N2 per unit O/H change is too large when compared with 
the data, for the pure secondary N models, while it is too 
small for the pure primary ones. 



The abundance parameter N2 vs. the ionization param- 
eter sensitive ratio (for ionizing temperatures higher than 
about 35000 K) log ([On] /[Oiii] ), is shown in Figure 2, 
together with the photoionization model results. The lines 
represent again photoionization models with the same ion- 
ization parameter and covering the whole range of metallic- 
ity. As can be seen, the lines for constant ionization param- 
eter are almost vertical and at a given metallicity the range 
in N2 covered by changing the ionization parameter is much 
smaller than the observed range. This confirms that most 
of the observed trend of N2 with O /H is due to metallicity 
changes. 



4 DISCUSSION 

Although relatively tight, the relation between N2 and abun- 
dance has a scatter that might be larger than the observa- 
tional errors. 

The comparison of models and data in Figure 2 sug- 
gests that some of the scatter in the N2 vs. 12+log(0/H) 
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Table 1. continued 



Galaxy 


12+log(0/H) 


N2 


Ref. 


Galaxy 


12+log(0/H) 


N2 


Ref. 


T1457-262 A 


8.170±0.074 


-1.362±0.065 


10 


NGC3310 E 


8.482±0.051 


-0.687±0.005 


21 


T1457-262 B 


8.230±0.081 


-1.666±0.058 


10 


NGC3310 L 


8.551±0.077 


-0.563±0.008 


21 


T1008-286 


8.174±0.073 


-1.471±0.078 


10 


NGC3310 M 


8.483±0.374 


-0.824±0.013 


21 


UCM1612+1308 


8.183±0.070 


-1.593±0.012 


15 


NGC7714 A 


8.473±0.051 


-0.509±0.011 


22 


C1409+120 


8.186±0.075 


-1.301±0.114 


10 


NGC7714 NllO 


8.553±0.179 


-0.305±0.044 


22 


N4212 C20 


8.192±0.090 


-1.093±0.027 


16 


NGC7714 B 


8.280±0.051 


-0.873±0.011 


22 


N4214 A6 


8.221±0.077 


-1.058±0.027 


16 


NGC7714 C 


8.265±0.167 


-0.859±0.070 


22 


N4214 C6 


8.374±0.087 


-1.197±0.027 


16 


NGC7714 N216 


8.607±0.199 


-0.322±0.036 


22 


T0440-381 


8.212±0.074 


-1.447±0.054 


10 


MIDI NGC5471 


8.182±0.119 


-1.640±0.157 


23 


Mrk 5 


8.214±0.090 


-1.309±0.053 


11 


MlOl NGC5471 


8.012±0.123 


-1.468±0.024 


19 


LMC 112 


8.225±0.084 


-1.612±0.061 


17 


MIDI NGC5471A 


7.986±0.142 


-1.553±0.054 


24 


Mrk 67 


8.226±0.094 


-1.721±0.068 


11 


MlOl NGC5471B 


8.295±0.146 


-1.095±0.019 


24 


T1116-325 


8.336±0.096 


-1.537±0.080 


10 


MlOl NGC5471C 


8.187±0.179 


-0.994±0.030 


24 


M33 CC93 


8.497±0.164 


-0.540±0.021 


18 


MlOl NGC5471D 


8.200±0.210 


-1.553±0.054 


24 


M33 IC142 


8.583±0.142 


-0.694±0.031 


18 


MlOl NGC5471E 


8.158±0.205 


-1.757±0.087 


24 


M33 NGC595 


8.407±0.049 


-0.816±0.012 


18 


MlOl NGC5455 


8.353±0.235 


-0.925±0.026 


24 


M33 NGC595 


8.461±0.188 


-0.735±0.029 


19 


MlOl NGC5455 


8.480±0.188 


-0.771±0.153 


23 


M33 MA2 


8.401±0.073 


-1.031±0.033 


18 


MlOl NGC5461 


8.499±0.339 


-0.854±0.022 


24 


M33 NGC604 


8.414±0.061 


-0.902±0.003 


18 


M51 CCM72 


9.095±0.215 


-0.493±0.013 


25 


M33 NGC604 


8.260±0.131 


-0.940±0.020 


19 


M51 CCM24 


9.040±0.222 


-0.464±0.063 


25 


M33 NGC588 


8.348±0.099 


-1.553±0.033 


18 


M51 CCMIO 


8.973±0.206 


-0.413±0.043 


25 


M33 NGC588 


8.395±0.244 


-1.346±0.023 


19 


NGC628 H13 


8.194±0.023 


-0.761±0.008 


26 


M33 IC131 


8.357±0.263 


-0.959±0.026 


19 


NGG628 H3 


8.234±0.024 


-0.784±0.024 


26 


NGC2403 VS35 


8.430±0.086 


-0.706±0.019 


20 


NGC628 H4 


8.310±0.026 


-0.679±0.040 


26 


NGC2403 Vb24 


8.423±0.129 


-0.828±0.019 


20 


NGC628 H5 


8.340±0.028 


-0.678±0.040 


26 


NGC2403 VS38 


8.365±0.150 


-0.860±0.019 


20 


NGC925 CDTl 


8.520±0.038 


-0.693±0.029 


26 


NGC2403 VS44 


8.332±0.125 


-0.876±0.019 


20 


NGC925 CDT2 


8.719±0.051 


-0.602±0.037 


26 


NGC2403 VS51 


8.353±0.196 


-1.012±0.020 


20 


NGC925 CDT3 


8.505±0.041 


-0.639±0.035 


26 


NGC2403 VS3 


8.333±0.192 


-0.948±0.019 


20 


NGC925 CDT4 


8.414±0.032 


-0.682±0.018 


26 


NGC2403 VS49 


8.362±0.157 


-1.106±0.021 


20 


NGC1232 CDTl 


8.440±0.032 


-0.441±0.012 


26 


NGC3310 Nuc 


8.831±0.059 


-0.298±0.007 


21 


NGC1232 CDT2 


8.597±0.041 


-0.602±0.024 


26 


NGC3310 A 


8.212±0.023 


-0.797±0.003 


21 


NGC1232 CDT3 


8.519±0.037 


-0.508±0.027 


26 


NGC3310 B 


8.441±0.032 


-0.638±0.007 


21 


NGC1232 CDT4 


8.516±0.037 


-0.578±0.023 


26 


NGC3310 C 


8.482±0.035 


-0.719±0.003 


21 


NGC1637 CDTl 


8.226±0.024 


-0.438±0.026 


26 



References to the table. 

(lO)Campbell, Terlevich & Melnick 1986; (ll)Garnett 1990; (15)Rego et al. 1998; (16)Kobulnicky & Skillman 
1996; (17)Matliis, Chu & Peterson 1985; (18)Vnchez et al. 1988; {19)Garnett 1989; (20)Garnett et al. 1997; 
(21)Pastoriza et al. 1993; (22)Gonzalez-Delgado et al. 1994; (23)Shields & Searle 1978; (24)Kennicutt & 
Garnett 1996; (25)Di'az et al. 1991; (26)Castellanos, Di'az & Terlevich 2001. 



relation could in principle depend on the degree of ioniza- 
tion of the nebula that in turn may depend on the age of 
the ionizing cluster. To cover the whole data range in the 
N2 vs. O/H relation, it was needed to vary the ionization 
parameter U from log(U)—-3 to -1.5, with a best fit corre- 
sponding to log(U) = -2. 

In an attempt to understand the dispersion in the O /H 
- N2 relation, we have studied the connection between the 
scatter of the N2 parameter and an age sensitive parame- 
ter, the equivalent width of H/3, EW(H/3). The ranking of 
EW(H/3) with age (Copetti, Pastoriza & Dottori, 1986; Lei- 
therer & Heckman 1995) goes in the sense that EW(H/3) 
in emission reaches values of several hundred A for an H ll 
region photoionized by a zero-age coeval stellar population 
and decreases as the cluster evolves and the stellar age in- 
creases. No significant trend was found between these two 
variables for the bulk of our sample, but a clear trend is 
visible for those with the highest quality EW(H/3) and O/H 
determinations. This result suggests that the presence of a 
range of ages in the data may be responsible for part of the 
scatter in the N2 vs. O/H diagram. This trend is shown in 



Figure 3 where we plotted the distance to the regression line 
in Figure 1 versus the EW(H/3). There is a clear trend sug- 
gesting that part of the scatter is correlated with EW(H/3). 
Intrinsic variations in the N/O abundance ratio will of course 
constitute another source of scatter. 

We can also expect dispersions of the order of 0.05-0.2 
dex in metallicity due to the differences between global spec- 
tra and smaller aperture exposures (Kobulnicky, Kennicutt 
& Pizagno 1999). This difference becomes more important 
when comparing low redshift galaxies with high-z spatially 
integrated spectra. Several effects my cause the emission-line 
ratios to produce oxygen abundance estimates significantly 
lower than those derived from small aperture observations 
of individual Hll regions. Kobulnicky et al. (1999) demon- 
strated through simple modeling, that temperature fluctua- 
tions are the primary cause for overestimating the electron 
temperature and underestimating the oxygen abundance, 
and that ionization parameter variations further exacerbate 
this systematic underestimation. We therefore consider our 
emission line measurements and abundance determinations 
from global spectra as lower limits. A more careful analysis 
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model log U = —3.0 



-2.5 -2 -1.5 -1 -0.5 

N2 

Figure 1. Oxygen abundance vs. the N2 calibrator (see text) for the whole dataset. The linear fit to the data, heavy solid line, taking into 
account the dispersion in the distribution, gives 12-|-log(0/H) = 9. 12-1-0. 73-N2, with a correlation coefficient of 0.85. The crossed-circle 
dots correspond to the photoionization model points for the metallicities as labelled by the side of the dots. The symbols refer to the 
galaxies in Castellanos, Diaz & Terlevich 2001 (CDTOl); six metal-poor galaxies from Terlevich et al. 2001 (Paper II), and a selected 
sample of low metallicity galaxies with references in Table 1. The only galaxies not listed in Table 1 come from the catalogue of Hn 
galaxies of Terlevich et al. 1991 (T91). 



of spatial variation effects on emission line intensities will be 
carried out in Paper II. 

The N2 parameter has clear observational advantages 
for ranking metallicities in starforming galaxies. Besides it 
being single- valued, it also appears to have a tighter corre- 
lation with O/H than R23; still, it has the drawback of using 
forbidden lines of N, making the abundance calibrator sen- 
sitive to variations on the N/O abundance ratio. Photoion- 
ization models taking into account primary plus secondary 
nitrogen encompass the data in the N2 vs. 12-flog(0/H) 
plane. To reproduce the dispersion it is enough to vary the 
ionization parameter U in the models by (-1-0.5,-1.0), with 
an average fit of log(U) = -2.0. 

We can expect that for the very low metallicity galaxies, 
the nitrogen lines would become less sensitive to metallicity 
due to the increasing presence of primary N in comparison to 
secondary N. Although we do not clearly see this behavior in 
our sample, we find that an alternative line ratio for the low 
metallicity branch would be useful, as the [N 11] lines also be- 
come too faint and difficult to deblend from Ha once in this 
regime. We adopted the ratio log([Sii] AA6717,673lA/Ha) 
(S2) as an alternative. In Figure 4 we show the relation of S2 



with metallicity. This relation still has the advantage of not 
being strongly dependent on reddening corrections, but its 
scatter is apparently larger than for the N2 - 12-|-log(0/H) 
relation. 

We emphasize that the N2 relation is clearly not bet- 
ter than the 823(4) methods, although it is observationally 
much easier to obtain, specially for moderate redshift ob- 
jects. It therefore has the advantage of allowing a very fast 
ranking of galaxies in a metallicity sequence, in particular 
when searching for low metallicity objects. A test for this 
method was performed on the Durham redshift survey data 
and is the subject of Paper II. 



5 CONCLUSIONS 

We have collected from the literature a representative sam- 
ple of spectroscopic measurements of starforming galaxies 
covering a wide range in metallicity (7.2 ^ 12-|-log(0/H) 
<, 9.1), and recalculated oxygen abundances in a self- 
consistent manner. We confirmed previously published re- 
sults on the correlation between the logarithmic ratio of 
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Figure 2. The N2 parameter plotted against the ionization sen- 
sitive parameter log ([On] /[Olll] ). The lines join models with 
the same ionization parameter for abundances [O/H] = 0.04, 0.2, 
1, 2 Zq. 



Figure 4. Oxygen abundance versus intensity ratio of the emis- 
sion lines [Sii] A6717A + [Sii] A673lA over Ha. The high metal- 
licity sample is from the compilations of Diaz & Perez-Montero 
(2000) and Castellanos, Diaz & Terlevich (2001); the low metal- 
licity sample is found in Table 1. 
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moderate resolution spectra, and the N2 line ratio does not 
depend on reddening corrections or flux calibration. 

On the negative side, N2 is sensitive to ionization and 
O/N variations implying that, strictly, it should be used 
mainly as an indicator of galaxy- wide abundances. 

Interestingly, the comparison with photoionization 
models indicates that the observed N2 is consistent with ni- 
trogen being a combination of both primary and secondary 
origin. 

Our main conclusion is that the combination of N2 and 
S2 provides a tool to roughly map the metallicity of galaxies 
(and search for low Z galaxies) using survey quality data, like 
e.g. that of the Sloan Digital Survey, for redshifts in excess 
of 0.2. 



-1 ' ' ' ' ' ' ' ' ' 

100 200 300 

EW(H|S) LAJ 

Figure 3. Residuals of the linear regression shown in Figure 1 
against equivalent width of H/3. The low metallicity sample can 
be found in Table 1. 



[Nil] /Ha (N2) and O/H and obtained an improved calibra- 
tor thanks to a larger sample of carefully and consistently 
calculated abundances and line ratios. 

Considering the weak points and limitations of empiri- 
cal abundance determination methods, we reckon that using 
N2 as a metallicity calibrator presents several advantages: it 
involves easily measurable lines that are available for a large 
redshift range (up to z ^^2.5); the N2 vs. metallicity relation 
is monotonic; the [N ll] and Ha lines can be separated in even 
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